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Abstract 

We investigate the electronic structures of some defective boron nitride nanotubes (BNNTs) un- 
der transverse electric fields within density-functional theory. (16,0) BNNTs with antisite, carbon 
substitution, single vacancy, and Stone- Wales 5775 defects are studied. Under transverse electric 
fields, the band gaps of the defective BNNTs are reduced, similar to the pristine ones. The energy 
levels of the defect states vary with the transverse electric field directions, due to the different 
electrostatic potential shift at the defect sites induced by the electric fields. Therefore, besides 
electronic structure and optical property engineering, the transverse electric field can be used to 
identify the defect positions in BNNTs. 
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I. INTRODUCTION 



Boron nitride nanotubes (BNNTs) are wide-gap semiconductors regardless of their diam- 
eter, chirality, or the number of tube walls.- This implies that BNNTs own uniform electric 
and optical properties,-'2ii^i^>^ and makes BNNTs suitable for nanoscale electronic and optical 
applications. The external electric field is useful to tune the electronic properties of nano 
materialSi^iS Khoo et al.^ found that, when transverse electric fields are applied, pristine 
BNNTs show uniform gap closure behavior. Transverse electric fields can move the con- 
duction band and valence band to two sides of the tubes. The band gaps of BNNTs are 
reduced almost linearly with the increase of electric fields. Under the same electric field, 
they are reduced with the increase of the diameter of BNNTs, independent of the chirality. 
The static dielectric response,- the screened polarizabilities,— and the optical properties of 
BNNTs under finite electric fields^i have also been investigated by theoretical calculations. 
On the other hand, the electronic structures and properties of BNNTs can be modified by 
introducing intrinsic defects.— i^^*^ The intrinsic defects can introduce localized states in the 
band gap.— Defects and impurities can be used to obtain n-type or p-type semiconducting 
BNNTs."!^ The role of defects and also the role of polarization field in BNNTs have been 
reported.— However, to the best of our knowledge, it is still unclear how transverse electric 
fields affect the electronic structures of BNNTs with intrinsic defects. 

In this work, we report the first-principles studies on the electronic structures of (16,0) 
BNNTs with antisite, carbon substitution, single vacancy, and Stone- Wales 5775 defects 
under transverse electric fields. Electric fields at three different directions perpendicular to 
the tube axis are applied. 

II. METHOD AND MODEL 

Our electronic structure calculations are performed by means of the spin polarized 
density-functional theory implemented in Vienna ah initio simulation package (VASP).— 
The projector augmented wave (PAW)^ method is used to describe the electron-ion inter- 
action, and the cut-off energy is set to 400.0 eV. Perdew-Wang functional^i is used for the 
generalized gradient approximation (GGA). The total energy is converged to 10^^ eV. The 
atomic structures are fully relaxed without electric field until the forces are less than 0.01 
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eV/A. The optimized geometrical structures are used to calculate the electronic structures 
under transverse static electric fields. Our test calculations on a (16,0) BNNT with B antisite 
indicate that geometry relaxation under electric field gives similar results. 

We adopt (16,0) BNNTs with diameter 12.8 A in our calculations. The tube axis is along 
the z direction. To minimize the interaction between adjacent images, we use a large cubic 
supercell (25.000 x 25.000 x 4.349 A^) for the pristine BNNT. For BNNTs with antisite 
and carbon substitution defects, we use a supercell with c=8.698 A, twice of the lattice 
parameter for the pristine BNNT. And for BNNTs with vacancies and Stone- Wales defects, 
we use a supercell with c=17.397 A. The special k-points for Brillouin Zone integrations are 
sampled using the Monkhorst-Pack scheme.-- For the above three models, 1x1x11,1x1x7, 
and 1x1x3 k-points are chosen, respectively. 

The method of applying static electric field implemented in VASP is in the spirit of the 
work of Neugebauer and Scheffler.— The electrostatic potential decreases along the direction 
of electric field. For convenience, all the defects are put at the positive x side in the supercell. 
The transverse electric field applied from the pristine side to the defective side is thus at 
the positive direction along x axis (S^). And the effects of electric fields applied at the 
positive y direction [Z^ and negative x direction (£_x) are also studied. The angle between 
the electric field and the positive x direction in the xy plane [ff] can be used to identify 
the direction of the electric field. The values of for £a;, £y, and are 0°, 90°, and 
180°, respectively. For (16,0) BNNTs, a large electric field is required to change the band 
structures significantly. The maximum strength of the applied electric fields in this study is 
0.3 V/A. The strength here is enough to induce distinct change of electronic properties of 
the tubes but no electron emission. It is expected that similar results should be obtained 
for larger BNNTs under weaker electric fields. 

III. RESULTS AND DISCUSSION 

First, we examine the electronic properties of a pristine (16,0)BNNT under electric field. 
The pristine (16,0) BNNT is a semiconductor with a 4.47 eV band gap. The conduction 
band mainly comes from B atoms, and the valence band mainly comes from N atoms. When 
the electric field increases from 0.0 to 0.1, 0.2, and 0.3 V/A, the band gap decreases from 
4.47 eV to 4.04, 3.52, and finally 2.96 eV. It changes almost linearly with the increase of 
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electric field strength, which is consistent with the work of Khoo et a/. The charge densities 
of the conduction band edge (CBE) states are moved along the direction of the electric field, 
to one side of the tube. And the valence band edge (VBE) states are moved to the opposite 
side of the nanotube. 

We consider seven kinds of defects on the (16,0) BNNTs: boron antisite (-Bat), nitrogen 
antisite (Nb), carbon substitution in a boron site (Cb) or a nitrogen site (Cat), boron vacancy 
(Vb), nitrogen vacancy (V/v), and Stone- Wales 5775 (SW). The results without electric field 
are presented first. The relaxed geometries are shown in Fig. [H The formation energies 
presented in Table [Dare calculated as 

Eform = Ebinding[def-tube] - Ehinding[ideal-tube], (1) 

where the Ehinding[ideal-tube] and Eunding[def-tuhe] are binding energies of systems with- 
out and with defects, respectively. The results are a little larger than those in the (8,0) 
BNNT.— i^i^i We also present work function of these BNNTs in the Table [H The band 
structures of the pristine and defective nanotubes are shown in Fig. [2], which are similar 
with the previous results. The defect A^^ introduces only occupied defect state, Vb intro- 
duces only unoccupied defect states, all the other defects introduce both unoccupied and 
occupied states in the band gap. In (16,0) BNNTs with and -Bjv, the unoccupied defect 
states are close to the VBE, but farther than that in the Vb case. While in the tube with 
Cfi, the occupied defect states are near to the CBE. The defect states in BNNTs with VJv 
and SW are deep in the band gap. The band structures of these systems can be sketched 
schematically in Fig. [3|, we define the energy difference between the CBE and the VBE as 
Ecv'i the one between the lowest unoccupied defect state and the VBE as Ep] and the one 
between the highest occupied defect state and the CBE as En. From Table [TTl we can see 
that, under zero electric field, the introduction of defects does not change the Ecv- As for 
the charge densities of the defect states, they are localized near the defect sites. It is shown 
from the profiles of the defect state charge densities in Fig. |H For BNNTs with carbon sub- 
stitution and single vacancy defects, there are also local magnetic moments in these systems 
(see Table UTTj) . 

For the defective (16,0) BNNTs, we study the effects of the 0.3 V/A electric fields at 
different directions. The Ecv in the defective systems are all reduced by the electric fields, 
with slight dependence on the field directions. In the narrowed band gap, the energy levels 
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of the defect states shift with the field directions remarkably, independent of whether the 
defect states are occupied or not. The 8,x lower the energy levels of the defect states. The 
8,_x uplift the energy of the defect states. The energy levels of the defect states are between 
the above two cases, when 8,y is applied. As shown in Table [TTl when electric fields are 
applied, with the increase of the angle (3, the relative movement of the defect states and 
the band edge states makes Ep increase and En decrease. In some defective systems, the 
lowest unoccupied defect state can be pushed into the conduction band by the electric fields, 
it makes the absence of the Ep in some cases in the Table HIl And the highest occupied 
defect state can be also pushed into the valence band. The localized charge densities of these 
defects keep almost unchanged under electric fields. In Fig. m as an example, the profiles 
of the charge density of the lowest unoccupied defect state of the BNNT with B]\f defect 
without and with 0.3 V/A E,-x are shown. 

Although the energy levels of defect states move almost linearly with the increase of P 
in all cases, it affects the concentration of carriers differently, depending on the kind of 
defects. For BNNTs with Vb and Cn defects, the 0.3 V/A S^, results in a small Ep and 
enhances the p-tjpe conductivity. The 0.3 V/A £„a, reduces the En and enhances the re- 
type conductivity most for the BNNT with Cb defect. In other situations, the Ep and 
En are maybe still too large to improve the conductivity significantly. However, we can 
use the change of conductivity under electric fields with different directions to identify the 
defect positions. If we change the directions of the transverse electric fields and measure 
the conductivity, for BNNTs with p-tjpe carriers, when the electric field is applied from the 
pristine side to the defective side, we would get the max conductance. For BNNTs with 
n-type carriers, we can get the max conductance at the reversed field direction. 

For the BNNTs with carbon substitution and single vacancy defects, the local magnetic 
moments show slight dependence on the field directions (see Table IIIip . it can be related 
to the change of Ep or En- Especially, in the BNNT with Vb defect, when the electric fields 
are applied, with the increase of (3, the magnetic moment is reduced. As shown in Fig. [2](b), 
in the band gap, there is one unoccupied defect state in the majority spin channel, and two 
in the minority spin channel. The reduction of Ep can increase the occupation probability 
on the lowest unoccupied defect state in the majority spin channel, so the total local spin 
can be increased, and the local magnetic moment can thus be enhanced. Other cases can be 
understood in a similar way. When the Ep and En are not so small, the magnetic moment 
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changes little with the electric field. 

To understand the behavior of defect states under electric fields, we present an explanation 
based on the field induced electrostatic potential shift. When £3, is applied, the defect is 
at the low potential side. While the field is reversed, the defect becomes at a relative high 
potential. So the electrostatic potential at the defect site can be related to (3. If we take 
the center of the BNNTs as origin, at which the electrostatic potential is set to zero, the 
potential induced by electric fields at defect sites can be expressed as 

AP = -|e| XXX cos{p), (2) 

in which |£| is the strength of the transverse electric field (0.3 V/A), x is the x coordinate 
of the defect site. As shown in Fig. [5l we examine the energy levels of the defect states and 
the band edge states with the AP induced by electric fields, taking BNNTs with and V/v 
defects as examples. The energy levels of the localized defect states are raised significantly 
with the increase of electrostatic potential at the defect sites. However, the levels of CBE 
and VBE are little affected compared to those of the defect states. The behavior of En and 
Ep under electric fields thus can be ascribed to the change of electrostatic potential at defect 
sites. And it can be expected that, stronger 8,x would bring down the electrostatic potential 
at the defect sites further, and reduce Ep more. The reversed stronger electric fields would 
reduce En similarly. 

In most situations, Eqv is reduced to the value similar to the pristine BNNT under the 
electric field. But in some cases, the Eqv is unexpectedly large. Those are BNNTs with Vb, 
Cjy, and B]y under 0.3 V/A £^ and BNNTs with Cb under 0.3 V/A £_^.. The large Ecv is 
related to the corresponding small Ep or En in these cases. The small Ep or En enhances the 
metallicity of that system, and thus increases the screening. With the stronger screening, 
the reduction of electrostatic potential along the direction of electric field is slower, therefore, 
Ecv reduction is smaller. It can be seen from the profiles of yz-plane-averaged electrostatic 
potential under 0.3 V/A £^. and 8,^^, as shown in Fig. [61 

From the above results, we can see that, the electric fields can adjust the electronic 
structures of the defective BNNTs in a different way from pristine BNNTs. The transverse 
electric fields narrow the band gaps, move the defect states near to the VBE or CBE, and 
can change the conductive and optical properties smoothly. Under the same strength of 
electric fields, the quantities of Ep and En change monotonously with the increase of the 
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angle j3. Compared to Ep and the Ecv changes httle with (5. So transverse electric 
fields can create shallow acceptor or donor states, depending on the kind of defects and the 
field directions. For one kind of defective BNNTs, one can measure the conductivity under 
different direction of electric fields, to identify the positions of the defects. The change of Ep 
and En under electric fields can also affect the local magnetic moments in the BNNTs with 
single vacancy and carbon substitution defects. The charge densities of the localized defect 
states are all almost unaffected by the electric fields. Besides changing their conductivity, the 
electric field engineering of the electronic structures of the defective BNNTs can also modify 
their optical properties. These external field tunnable properties facilitate application of 
defective BNNTs as nano devices. 

IV. CONCLUSION 

In summary, we study the electronic structures of defective BNNTs under transverse 
electric fields. The band gaps of these BNNTs can be reduced by electric field like that 
in pristine BNNTs. The relative positions of the defect states in the band gap are moved 
when the transverse electric field direction is changed. The charge densities of localized 

defect states are not affected by the electric field. The behavior of electronic structures 
can be ascribed to the difference of the electrostatic potential induced by electric fields at 
the defect sites. Our results indicate that defective BNNTs can be used as electric field 
controlled nano electronic or optical devices. 
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TABLE I: The formation energies {Eform) and work function ($) of seven kinds of defects on the 
wall of the (16,0) BN nanotubes. 

Bn Nb Cb Cn Vb Vn SW 
Eform{eV) 5.37 7.34 2.70 1.55 16.08 12.01 5.71 
^{eV) 4.86 3.91 2.58 5.33 5.74 3.76 4.81 



TABLE IL The values of -EcVi and En as plotted in Fig. [3] of pristine and seven kinds of 
defective (16,0) BNNTs under zero and 0.3 V/A electric fields. The electric field is applied at 
positive X, y, and negative x direction, respectively. 
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TABLE IIL The magnetic moments per ceh (M) of (16,0) BNNTs with Cb, Vb, Cn, and Vat 
defects under zero and 0.3 V/A electric fields. 

M{hb) £^ e.y 

Cb 1.00 1.00 0.99 0.90 

Vb 1.04 1.41 1.27 1.02 

Cn 0.99 0.87 0.97 0.99 

Vn 1.00 1.00 1.00 1.00 
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FIG. 1: The optimized structure of BNNTs with defects. Pink ball is boron, blue ball is nitrogen, 
and carbon is in gray. 
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FIG. 2: Band structures of BNNTs with defects (a) A^b, (b) Vb, (c) Cat, (d) Bjq, (e) C^, (f) Vat, 
(g) SW and (h) pristine BNNTs without electric field. The red is the majority spin channel, and 
green is the minority spin channel. The dash line represents the fermi energy. 
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FIG. 3: The schematic diagram of the energy gaps between the defect states and the band edge 
states. Ecv is the gap between the conduction band and valence band edge states; Ep is the gap 
between the lowest unoccupied defect state and the valence band edge state; and En is the gap 
between the highest occupied defect state and the conduction band edge state. Ef is the fermi 
energy. 
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FIG. 4: The profiles of the F point charge densities of the highest occupied defect states of BNNTs 
with (a) A''^, (b) Cb, and (c) Vn', and the lowest unoccupied defect states of BNNTs with (d) Vb, 
(e) Cn, and (f) SW under zero electric field. The axial view and side view of the lowest unoccupied 
defect state of BNNT with Bn (g) under zero electric field and (h) 0.3 V/A £.^x are also presented. 
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FIG. 5: The energy levels of the valence band edge (VBE) state, the occupied or unoccupied defect 
states, and the conduction band edge (CBE) state of BNNTs with (a) Bj\f and (b) defects. 
The abscissa is the the shift of electrostatic potential at the defect sites induced by electric fields, 
calculated with Eq. [2j The strength of applied electric fields is 0.3 V/A. 
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FIG. 6: The profiles of yz-plane-averaged electrostatic potentials of pristine and seven defective 
BNNTs under (a) £^ and (b) £_^. 
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